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FOREWORD 


This report covers a specific section of the Statement of 
Work, Exhibit "C", of Contract NAS9-9200 concerning tide-produced 
gravity gradients and crustal stresses on the moon and their 
possible link with lunar transient events and moonquakes. This 
is one of four different areas of investigation performed under 
Exhibit "C" and the results have been prepared as a separate Part 
II of the Pinal Report for convenience of future reference to the 
different subject areas. This Part II and the three reports on 
the other work areas (Parts I, III and IV) constitute the Pinal 
Report due under the Contract, 

The results required by the Statement of Work, Exhibit "C", 
are essentially self-evident from the mathematics, calculations 
and analyses performed and presented in this report. All items 
of the Statement of Work are believed to be more than adequately 
satisfied. 

The role of the tide-produced gravity gradients as triggering 
mechanisms for moonquakes and lunar transient events has long been 
a subject of interest of the junior author, Lloyd G.D. Thompson 
and the Contract Technical Monitor, William B. Chapman. The 
preliminary investigation covered by this report was chiefly 
performed by Mark H. Houston under the direction and guidance of 
Thompson and Chapman. In recognition of this, Houston has been 
selected as the senior author. 
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PART II 


LUNAR TIDAL GRAVITY GRADIENTS AND STRESSES 


AB, .T 


Preliminary analysis of the gravity gradients associated 
with gravity tides on the moon caused by the earth indicates that 
the relative changes in the gradients are very Irregular and 
large and about 15 times greater than those experienced on earth. 
Thus gradients, in preference to gravity tides themselves, may 
well be an important key in correlating tide effects with lunar 
transient events and moonquakes and also in determining trigger- 
ing mechanisms for crustal movement and faulting. Preliminary 
analysis of lunar crustal stresses and straint.^ caused by lunar 
gravity tides indicates that these factors may be more direct 
causative agents or triggering mechanisms. In particular, the 
cubic dilation undergoes relatively large changes and is about 
11 times greater on the moon than on earth. Thus it should be 
correspondingly more Important, Development of formulae for the 
gravity gradient tensor and the stress tensor terms plus computer 
programs for calculating lunar tidal gravity, gradients, stresses 
and strains together with some suggested ways in which these terms 
may play a role in crustal mechanics provides a starting point for 
further more detailed investigations. 
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T. Tn1:rodunV, 1 on 


r.unar tr/m;’ loni ovontn ocourlm’; a.i f*.'d patclKo'. , tirlphi np'oin: 

In iiiiadov/od aroan ami other phonomena have been ohiujrvod on the .’.ur- 
faoo of the moon for over three liundred yearii . Vlhethor then.e ev<mt:’ 
are volcanlfi eruptiom.i, paaeoun outflows, id.atio-aharf^ed du;;t cloud;’, 
nr merely optical effocta la atill unknown, 

dovoral Invontlpatnra have aaaerted that trannient lunar evmit;’. 
may be trimRoi-’od by tidal otrei’n^ec within the lunar cruot and have 
examined the correlation between the variation of tiie vert i cal f^ravity 
tide and the incidence of lunar ovent;">. Thl;’ report developer and 
presents the formulae for the f?;ravity p;radlent c and tidal otrenr.en 
accociated with the gravity tides on tlie moon caused by the earth. 

It further indicates uhe role of these term.s an possible trij’if’ierlnn 
mochanislms for lunar transient events or moonquaken , This work is 
intended to form the basin for more detailed Investigations to estab- 
lish a correlation between tidal gravity gradients and stresses and 
lunar events and the role of these terms as triggering agents . The 
equations and analysis have been considered for the case of gravity 
tides on the moon caused by the earth only with the solar effects 
omraitted. Basically, the formulae also apply for the gravity tides 
on the earth caused by the moon but in this case the ephemerls of the 
moon must be substituted. 

First the gravity gradient tensor l.s developed and the various 
terms examined. Assuming particular forms and values for the Love 
Numbers A and h and the material constants 1 and y, a method for com- 
puting the crustal stresses using crustal strains derived from the 
tidal potential is presented. The tidal stresses and cubic dilation 
are then examined. 
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i J . TjUnar Tidal (iravlty Qrad l.ontf. 

a) 'ride Potential 

Jt can be ahown that in selenot’;i’aphio apiierlcal coor’dinattia 
the gravitational potential, U, on tlje surface of i' laoon can 
be expressed as: 

0 = ^ ^ (. 0 .. » _ 

where y ~ universal gravitational constant 
= mass of the earth 
d = earth-juoon separation 
I' = radius of the moon 

0 = included angle between the suh-terrestial point and 
the observation location 

In astronomical coordinates the included angle, 0, iuay be related 
to selenographic coordinates by 

cos 0 = sln<{) sin5 + gqs6 cos^ cqs^ cos! + sin^ sinl 

where 

(j) = latitude of the selenographic observation point 
6 = latitude of the eai’th (declination) 

X = longitude of the earth (hour angle) 

^ = longitude of the selenographic observation point 
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1'MiKiii.y , ..ubi'.tltutinr: the ex{Ji‘<?;in Lon toi' 0 Inicj rhie tlUa.i 
potential anh r'eaia'an;':lih;: Into ijaplacie',.; ;;epai'at l,on , ihio tidal 
potential may bo vn-itton 
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b ) gravity Tldea 

liari’ldon (Xyb3) and Sutton et al (1963) have dlacuoiied 
aapeuts of the tidal vai'latloiu; on the inoon. 

The lutiar gravity gradienta are beat considered In local 
coordinates In the lunar surface. If we adopt axes with u posi- 
tive upward^ x positive eastward and y positive nor'thwai'd, idien 
tire gravity forces in the Xj y, and a direction will be pjlveri 
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c ) Gravity Gradient Tensor 

similarly, tho components oT the gravity gi'adlenb tonsoi’ 
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A oornpnt'.'T’ ppoj’-vam for oaTnul at ! nr: tlfia'l. /’’r'nvlty ami iri'-'ivlty 
Nl’afllont (toinponnuti’ on i:h<» iriooh <luo to t)j" oaiM.li t;; f’llyr-ii hi 
Amu'iidl X A, 

Tn ronoral, f.h<? rradlont.;’ on l/no moon vary In a coitipli*x 
manner rivor tlio iiurraro of tin- iwm* bo<ly and tlm vari al; Inn;’ ar'.‘ 
not an rerinlar an timno on oartli. ':i’hn lon/tltudn of tJio “ar'th In 
tho lunar nky ohanfon nlov/ly wttb lunar 1 Miration and tlK? latltudo 
ohanfnn; with tho rovolutlon of tho moon In Itu orliit;, Additional 
variation In tdio uradLontr. In Introdurc'd by tho ulmn/>]l nil oartli-moon 
dir, tan GO . 

Ti:o variation in tho variou;! r.radlont componontu for' a lunar 
month period, tho GorriTinnont dlfforenGon one to another and tho dlff- 
or''"jncon from place to plaoo on tho moon are Illustrated in Nh’:ure;; 

1, f, and 3. The tidal contrlbutlfui to tho various f^radlont terms, 
ran bo soon In these f iaiures . Flp;ure 1 Illustrates the oharaator 
of tlio r;radlents of the vortical oomponent of pravlty for a hyrio- 
thf'tloal oa;'.e In which iho r;radlents. are obsorvod at a planetary 
latitude of and the t Ide-producinn body Hes. at a latltud':? of 

20'’s. For comparison, Flp.ure ,? illustrates, a simll.ar ca;m* In which 
tho latitude of the disturbini’; body has been clran/.iod to Lh/'^K, In 
addition, Fif’iure 3 shows an oxamnlo of the calculated /’iradlents 
11 „ d,,,, and tl for i.he crater Aristarchus. 

The vortical praadlont comrionent acts in the local vertical 

« t <-* 

and n,an be seen in Flruros 1, ;? and 3 to shift in r.'hase anple and 
magnitude at different sites on the lunar .rurfaoo. and ar^: 
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On tdo moon t,h(> maximum noak-*iand>'.*ai: variation In idr; voud, 1 o.a 1. 
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about 0.3 mral and the aorronpnnd ino nlianpio in bho /'I'ad i'.'nt In aln>ub 
b X r;al/fim or O.OOOh Eotvou Un1bi'. 


tt 1. I'lotvo.u Unit = I K. U. = I X 10 


-n 


f’;al/om = 1 


10 


-0 


cm/i 


' o c 


/am , 


7 



Local Angle of Earth (Degrees 

Variation of Grradients of Vertical Icmoor.e: 
Lunar Tidal Gravity at 4£'’D, Earth at IfS 





)f Gradients of Vertical Conponent of 
. Gravity at Earth at 20 





IXI. 


Lunar Tidal Stresses 


While the variations of the lunar gravity components and 
the components of the gravity gradient tensor may yield some 
correlation with the distribution of lunar transient events ^ 
further Insight into the problem may be obtained by considering 
tidal stresses induced within the lunar crust to be the trig- 
,.,ering mechanism for the lunar events, uenee it is necessary 
to model the response of a planetary body to a varying external 
gravity field. 

Dix (1968) j using Takeuchi's (1950) model of the earthj 
attempted to show that tidal stresses within the earth may be 
sufficient to provide a triggering mechanism for the earthquake 
in Kern County, California on July 21, 1952. In our analysis of 
stress in the lunar crust we will use a much simpler model for the 
moon, uur lack of knowledge of the interior structure of the 
moon precludes any complex model. 

a) Displacements 

Adopting the previously defined coordinate system, the dis- 
placement at the moon’s surface can be written, 

^x g cos^ 3X 

q - i M 

y s 
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v^liore 



U = i’^ravitational potential 
ii; = acceleration at the moon'.i surface 
h,£, = Love Number's for the vertical and horisontal tides 


respectively . 


Lubstltutinc the expression for the potential and asamiiln;;: 
g varies only with r gives j 
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b ) Strain Tensor 

To calculate the strain tensor, we will make the assumption 
that hj % and g are functions only of the radius r (s direction) 
and tnen proceed as before. The elements of the strain tensor are, 


e 


X 


'XX r COS0 dX 


1 

'yy " r 30 


S S, 

^ tan0 + ^ 


12 




3-. 


Q 

55Z 

_ /t 

3r 



_ 1 

+ ■ 

yx 

r 3(|) 


1 

as 

a 

®xz 

r costj) 

ax 


s: — ..jI, ^ 

3r 

fy. 

r 


, as 

• + J- n 

’x r cos(|) ax 


33 

X g 

^ r 3 (|) 


e = e 
xy yx 


e = e 

fp v V ^ 


e = e 

yg gy 
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We can reduce the expressions further by as Burning that the 
density and deformations are homothetlc with regard to the 
center. We assume that 
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G ) Mtreiiu Tcnsur 

Given tho above atrain teisoor ami aasumlin;; Hooke 'a Law 

for’ a liornogoneoua iootropic body, ahe otreaa tenaor* oau be 
written in generalised form an. 


^ '^IJ ®11 ^ ^'ij 


where X and y ai’e Lame' constants commonly known as the com- 
pressibility and rigidity. I’inally, since the general coordinate 
transformation of a second order tensor is written. 



^ia '^’ij 


the stress within the lunar crust along any direction, n, v/hose 
direction cosines are 2.^, and 2,^ are given by 
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! V . OooT^liy>' I ' ml.' m. - for IiUna i* l';V«*n 

f'ovMi'al pnn;vlblf cau;;<'ii'.1 vo a;-^ontt’. for lunar traruilont ijvent;; 
may bo oonn 1 ilorod . Unloua iiomo inochanliim oxlatii whloii acoumulatoi’ 
biio ildal onorry nvor lone porlo<ln of time, tho tidal onorey involvod 
i.5 much too ninall to lio a diroct cauuatlvo aermt . If faulting or 
moon quakor aro ronporurlble foi* the ob:’orvod phonomona, tiien tidal 
rtroaoLKi may tiervo ats a trleeorlne aeont , The provioua analy;'!;' 
ya’cviden a moano to calculate tidal utreisseL] In any direction and 
will thua allow an analyai'a of the tidal intreoM alone a particularly 
oriented fault plane, A particular method for calculating the 
tidal Btrecc alone <o fault plane is given in Appiendlx B, 

If the lunar crust contains reservoirs of inoompress ii)le fluid 
such as water or molten lava, or contains structures of widely 
varying compressibility, then variations In the cubic dilation of 
the lunar crust could be responsible for the transient events . On 
earth fluid levels in oil and viater wells are observed to change witli 
the phase of the solid earth tide. The volume of the liquid remains 
constant but the volume of the enclosing cavity varies with tlie cubic 
dilation of the crust. The tidal cubic dilation can be written 
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ouppoae now we examine the I'eiative ma/hiituUe of aiinpie 
extension and of cubic dilation on the moon as coirif'ui'ed to tiiat on tie 
earth, in order to do tills, we must select suitable hove Numbers 
for the moon and the earth. The cubic dilation, ii, can be 
written in terms of the potential ami a charactei'istlc number, f. 


u 



where f = r + bh - bH 
dr 

Similarly, crustal extension in the 'S direction is, 

ijU 


The two constants for the earth have been experimentally 
determined to be O.620 ( f „ ) and O.SB^l (h_) {Melchior, iybb, 

S 0 

p, 300}. Melchior (iy66, e, 32‘j+) has shown that for a homopeneous 
self gravitating sphere. 
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ID 


li 


ttie moon, 

" VjO K,aia 

1 ' = I.Y;-iH X 10° cm 

ami if we aasume that the rigidity of the moon 1 l 5 equal to that 
of tile earth 'j oruatj 

a = 10^^ {cgfi) 

In that caae, n = .Obb, t = .00/ and f = .0b!3. 

Now taking the ratio of the moon'.; cubic dilation to that 
of the eartii. 






Thu;? v/o !ic?e tliai wh U.o Uio radial nxtoii;? ! <>n on I ho moon 1;? only 
1.‘3 t?linoi? that; of tho ourth, tho amp 1:1 tudo of lunar ouhlo oxtjan;;Ion 
:1 ;? elovon tljiion that, found on oarth and .?hould In' aorrin?pt)nd ! Uf'.l ,v 
mnro Import.ant In lunar nrui'.tal moahanlan. 


To furtlior 1 iivo;?t 1 jiato nnd IllUn tratt? tht;? ormoopt , tlio monthly 

var’lat;lon in tho vortical and horizontal i?tro.n?(,';? and cubic dilation 

on tho moon duo to tho oarth have boon oalnulalc'd for tlio ;?amo tlino 

at t,he throo orator;? of Arl;?tai’Ohu;? , I'latn and Alf.iiioiu'.u;? arid aro 

;?hown in Flf>iui’o£; , ’j and 6 ro;?peotl voly , Tho otjinputor program for* 

porformlnf’: tho;?o caloulat:lon£? .1;? Rivon ,in Appondlx C. P.l iparf;;; J| , h , 

and A ahow that tho variation in tho cuhlo dilation ha;? a vo^ry larf:'» 

amplltudo while tlio vortical (T ) and liorizontal (T ) c,trr?;?;’,o;; 

XX 

exnorienco much nmallnr chanp;o;? . Furthor, the behavior of tho.'o 
quantitio;; i;? clf’inlf Icantly different at different location;.? on tin? 
moon. At Ari;.;tarohu;? and Pl.ato (Plfmrer? tl and h) tlr? cubic dilation 


and T have much the ;:.amc behavior but T„„ in ponitive at Ar l..?tarcliu 

Y-x /i XX 

and ncRatlve at Plato. At Alnhonnu;? (Fir;uro 6) the amplitude of the 
cubic dilation in con:?idorably reduced hut tho marnitude 1.? at a 
very hifih pos'.itive level. In addition, both T and T , are al;'o 
positive. ThlL? indicates that the 3trcn?;?e;? and particularly tlio 
cubic dilation could well be important cau;?atlve arentc in lunar 
cructal meohanict? and fault Inr;. 
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Figure Stress and Cubic Dilation 


at Crater Aristarchus 



Time (Days ) 


Figure 5* Stress and Cubic Dilation 
at Crater Plato 
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Figure 6. Stress and Cubic Dilation 


at Crater Alphonsus 


V. ConclualoiiH and Hecommendations 


Prom other work, there is an apparent correlation between 
lunar transient events and/or moonquakes and gravity tides which, 
by symmetry of equations, applies also to the tidal gravity 
gradients, Tidal variations in gravity are not likely, in them- 
selves, to be direct causative agents since the tidal energy 
involved is too small. On the other hand the relative changes in 
the gradients are significantly Irregular and large and are, in 
fact, 15 times greater on the moon than on the earth. This 
indicates that gradients (one or more components) could possibly 
relate more directly to lunar events and the generation of large 
local differential forces and therefore could be a significant 
factor or key in crustal mechanics. They may contribute to or 
have a modifying influence operating in conjunction with other 
physical mechanisms. More particularly, there may be a mechanism 
for accumulating tidal effects over a long period of time or there 
may be large forces created differentially when gradients are 
considered over large areas or across a fault zone. This prelim- 
inary work provides the basis for further and more detailed in- 
vestigations of these and other approaches to the problem. Such 
additional studies are highly recommended. 

Lunar crustal stresses and cubic dilation caused by gravity 
tides and gradients could well be important and more direct 
triggering mechanisms or causative agents in crustal mechanics and 
faulting. In particular, the lunar cubic dilation exhibits large 
variations and is eleven times that experienced on earth. Thus 



It should be correspondingly more important in lunar crustal 
mechanics. Further investigations of long term stress build-up 
and the stresses and cubic dilation generated over large adjacent 
areas and across a fault sone are logical and anticipated con- 
tinuations of this preliminary work. Such studies are highly 
recommended, not only In support of lunar Investigations but also 
in the interest of earthquake prediction to which this work also 
applies . 
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APPENDIX A 


Computer Program for Calculating 
Tidal Gravity and Gradient g on the 
Moon due to the Earth . 



n n o 


C 
C 

c 
c 
c 

IMPLICIT DOUBLE PRECISION (A..H»L,0-Y) 
dimension JYSAR{2/26) 12) 
data JIEAK{0)/2«l5019/» ( JDAY( 1 ,K) ,Kr:i , 12)/ 
10,31,60,91, 121, 152, 18 2,2 13, 2« 4, 274, 305, 335/ 
COSQ(THETA) = DCOS(THETA)*DCOS(THETa) 
SINQ(THETA) a DSIN(THETA) *DSIN (THETA) 
C0S2(THETA) = DCOS(2.0D0*THETA) 
SXN2(THETA) = DS1N(2.0D0*THETA ) 


LUNAR TIDES AND GRADIENTS 
M HOUSTON 5/14/70 




this program calculates the tidal anomalies and tidal 
gradients on the moon due to the earth, 


EMASS = 5,975U0 
GRAVCN = 6.67U0 
EMDIST = 3.84411D0 
RMOON s (2.73D-1 )t.(6.37l23D8) 

CONST = ( 1 .5D0) *GRAVCN*EHASS/(EMDIST*EMDIST*BHDIST) 

CONST a CONST* ( 1 .0D-1 1 ) 

JULIAN calendar 4/2/70 

DO 1 J=1,26 

1 JYEAR(J) a JYEAR(J-.1) + 1461 
DO 2 Ia2,4 

DO 2 Ja1 , 12 

2 JDAI(I,J) a JUAI(I-.1,J) + 365 
JDAY('2,1) a JUAY(2,1) + 1 
JDAI(2,2) a JDAY{2,2) + 1 
TYPE 99 

99 F0RHAT(/' number OF CASES (12)'/) 

NCASE = 0 

accept 104, ICASE 
999 type 100 

100 FORMaT(/' day, month, yeah ^ 312'/) 

WCASE a NCASE + 1 

accept 101,IDAY,IMONTH,1YBAR 

101 F0RMAT(3I2) 

11 = IYBAR/4 

12 = -4fI1+IVEAR+1 

JULdaY = JYEaK(II) + JDAY(I2,IMoNTH) + IDAY 
type 102,JULDAY 

102 FORHaT(/' JULXAN calendar day ai,Il5) 

type 103 

103 FORMAT(/i number OF DAYS TIDES AND GRADIENTS DESIRED = '/) 

accept 104,NDAY 

104 F0RMAT{I2) 

NDAY a NDAY+JULDAY 
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9id 

91 

92 

2«0 

201 

202 

190 

69 

C 

C 


TYPE 90 

F0RNAT(//' lunar observation POINT'//' LATITUDE s '/) 
accept 91,0BLAT 
FORHAT(U) 
type 92 

FORMAT!/' LONOITUCE « '/) 
accept 91fOBLONG 

type 200 

FORMAT(/i output file unit designation (12} = '/) 
accept 1014,IUNIT 

type 201 

FORMAT!/* OUTPUT FILE NAME (A5) = '/) 

accept 202,ZNAMB 

F0RHAT{A5) 

CALL OFILE!IUNIT,ZNaHE) 

WRITE (lUNIT, 190)OBLAT, OBLONG 

FORMAT!///, 1XJSHLUNAR OBSERVATION POINT COORDINATES,// 
1 7X8HLAT1TUDE, / X9 HLONGITUDE , // , 5X 1 PD 1 3 , 6 , 2X , D 1 3 . 6 ) 


RDEG = 1 ,7tl5329251994329D«2 

oblat >= oblat + hdeg 

OBLONG = OBLONQ^BDEG 
JULDAY = julday ~ 1 
julday = julday + 1 


lunar EPHEMEHIS U/20/70 


T = DBLE(FLOAT{ JULDAY-2415019) )/3.6525D4 

PI = 3,14159245358979300 

AI = !9.21D1)*PI/!(1.aD2)<f{6,0D1)) 

S = 2.704365902 + ( ( ( 2 . 0D«6 ) ♦T+ 1 , 9 8D-3 ) #T+4 , B 1 267 8 9057DS ) *T 
P = 3.343295602 + { ( ( -1 . 0D-5 ) *T- 1 . 032D-2 ) *T+4 . 06903403D3 1 *T 
H = 2.796966802 •*, ( ( 3 , PD-4 ) *T + 3 . 600076 B92D 4 ) #T 
RMIN = 2.9088U20866572l5Dr,4 
S - S*RDBG 
P = P*RDBG 
H = H*RDE0 


DRATIO s 1.0D0 + (5,45D~2)'<‘DC0S(S«P) + ( 2 . 97D-3 ) ^DCQS ! { 2 . 0D0 ) * 
1(S-P)) + ( 1,002O^2)*DCOS(S-Pr.2,0D0+(H-p) ) 

2 + (8.25D-'3)-hDCOS(2.0D0#(S-H) ) 

PS = 2.8122083D2 + ( ( ( 3 , 0D- 6 ) *T + 4 , 5d- 4 ) *T+ 1 , 7 1 90 2D0 ) * T 
pS = PS*RDEG 
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203 

204 
C 


Si - (3.77D2) *DSIN(S-P) + ( 1 , 3D 1) *DSIN ( ( 2 1 0D0 ) * { Sr,p ) ) 

1+ (7,6D1)*DS1H(S«P-(2,0D0)o(H-P) ) ~ { 1 . 1 D 1 ) *DS1N { H-PS ) 

2+ (4,0d 1 ) *DiiIN( (2,0D0)*{S-H) } 

Si = S + S1*Hl1JN 

an = 2,591832402 + { ( ( 2 , 0D«6 ) *T + 2 . 08D-3 ) 1 , 9 34 1 4 20 1 D3 ) *T 
an s AN*BDEG 

A = DSlN(AI)>fDC0S{S1-AN) 

AMD e (1,03P-2)-fDSIN((2.0D0)*(51-ANn*HDEG 

beta = (3.08n2) *DS1N(S-,AN) + ( 1 .7D1 ) ■«DSIN(2,0D0*S-P-AN) 

1+ ( 1,0Dl)*DbIN( (2,0D0)*(8-H)-S + AN) - ( 1 , 7 d 1 ) *DSlN ( P-AN ) 

beta = BETAtiHniN 

TANB = -DSTN ( AI) *DS1N (S1«AW)/DC0S (AI) 

ELAT = (DATAN(TANB)-bETA) 

ELONG “ (S1-S + AWlf + A*ELAT'/ 

TANB = -DSIW {AI)*DSIK(S1-AN)/DC0S(AI) 

ELAT “ (DATAN; TANB) -BETA) 

ELONG = (SI-s + AMU + Am-ELAT) 

E1 = ELAT/RDEG 

F2 = elong/rdeg 

WHITE (IUNIT,203)JULDAY 

FOBMaT(/////1X46HEARTH SELBNOGRaPKIC coordinates for JULIAN DAT^ 
1 110 ) 

WRITE (lUNIT, 204 )E1 ,E2,DHATI0 

F0EMAT(//, 7X8HLATlTUDE,7X9KLONGITUDE,6XaHPARALLAX, 

1 //,5X» 1PD13.6^2(2X^D13,6) ) 

TIDE AND gradient CALCULATIONS 4/28/70 

DRATIO sr dratio*dratio*dratio 

P = OBLAT 
D = SLAT 

L “ ELONG " OBLONG 

UZ = C0SQ(P)*C0S0{D)*C0S2(L)+SlN2(P)t'SlN2(D)*DC0S(L) 

1+ ( (3.0D0) kSINOCP)-! ,0D0) *(S1NQ{D)-1 .0D0/3 ,0D0) 

UZ s -rmoon»dratio»uz*const 

UY = -SIN2(P)*COSQ(D)*COS2<L)+COS2(P)*SlN2(D)*DCOS(L)*t2,0D0) 

1+ SIN2(P)*( t3.0D0) *S1N0(D)«1 ,0D0) 

UY = -HMOON*DHATia'«UY*(5,0D-1 )*CONST 

UX = DCOStP)*COSO{D)t‘SlN2(L)+DSlN(P)*SlN2(D)-fDSIN{L) 
ux = RM 00 N*DRATIOhUX*C 0 HST 
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U2Z =; UZ/HMOOW 


300 

301 

70 

71 

302 


=; SINg(P)*COSQ(D)*COS^{L)-SrN2(P)-<SIN2(D)*DCOS(I.) 

1+ < 3.0D0*COSQ(V)~1 ,0D0)*(SINO(D)-1.0P0/3.0D0) 
yYY = -DRftTlO*COHST*UltY 

UXX s COSU(tl't'COS2{L)+SlNg(D)-1,0D0/3,0D0 
UXX s DHATI0«C0NRT#UXX 

UYX = -DSIN(P)^C0SQ.(D)*sIN2(L)+DC0S{P)*SIN2(D)'»DSIH(L) 

UYX S (2.0D0) »DRATI0*C0KST^UYX 

UXZ = DCOS(P) *COSQ(D) *SXN2(L) -COSQ(P) *COSg(D) *C0S2(L) 

1+ SXN2(P) *SXN2(U) *(DSIN(L)"DC0S(L) ) -DSIN ( P } *SIK2 ( D ) *DSIN U) 

2- (3.0D0*SIKQ{P)-1 ,0D0)*{SIHQ(D)"1 ,0D0/3,0D0) 

UXZ = DBATIOi^CONST*UXZ 

UZY s -SIH2 (P)«kcOSQ(D}*COS2(L)+2,0D0#COS2{P)*SIN2(D)*DCOS(M 
1+ SIN2(P)*(3.0D0*SINQ(D)-.1.0D0) 

IIZY = -DRATI0*C0MST*1J2Y 

WRITE (IURIT,300)UX,l)Y,UZ 

P0HHAT(///, 1X30HTIDAL GRAVITY ADOHAllES (GALS),//, 

1 10X2HUX, 13X2HUY, 13X2HUZ,/»5X, 1 PD 1 3 , b , 2 ( 2X , D1 3 . 6 ) ) 

WRITE (IUNIT = 301,1 UXX.,UYX..UYY..tJXZ,,UZY,,UZZ 
POHMAT(///, 1X29HTIDAL GRADIENT MATRIX (SEC-2)^//, 

1 10X3HUXX,/,5x1PD13,6,//, 10X3HUXY, 1 2X3HUY Y , / , 5X 

2 2(D13.6,2X) ,//, 10X3HUXZ, 12X3HUYZ, 1 2X3HUZZ , / ^ 5X , 3 { D 1 3 , 6 ^ 2X ) ) 

IP (JULDAY-NDAY) 69,70.70 
CONTINUE 

IF (NCASB-ICASE) 999,71,71 

continue 

type 302 

FORHaK//' FINIS') 


STOP 

pND 
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APPENDIX B 

Mothod for Calculation of Tidal dti’oc.i3e 
Alonp, 0 , Fault I’lane . 


3f 



CALCULATION OF TIDAL STHlynS ALONG 
A NAULT PLANK 


The following is a method for resolving the tidal stresses 
along a fault plane whose normal is denoted by tiie Uireetion co- 
sines and x^). Tlie formulation transforms the 

stress tensor into a new coordinate system which contains 
e^j ' . The stress normal to the plane and the aalmuthal 

dependence of the stress within the plane oan then be easily 
computed. 


1) Normal to desire plane of calculation has direction co- 
sines 

2 ) Tilt new coordinate system so that x., ' (a) has direc- 

tional cosines of l.o. 



A 


1 




A 


3 


3) Rotate new axes around z' axis until the angle between 
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x' unii ;• ;ixir, 13 it/;' oi* 



,V( 




For the transformation of a second order tensor e, , the 
transformation is^ 



^nj ^inn 

3h 


I 



'i'hu tonal on in any dii’octloji (with 

dinutoa) v/itli dlrootion cooines a-,, a,,* 

X d 


I’oapect to the new ooor- 
will be fii^iven by 




“ii 


a. 


iiti 


+ a. 


■".SZ + “l“2 


* “l“3 ■'‘13 “2“3 ‘•‘23 


at the aurface of the moon because of limiting conditions 

Along the normal to the defined tilt plane which is ' in 
our new coordinate system^ the stress will be, 


Finally, stress in the tangential plane of the surface of the moon 
as a function of azimuth will be, 

= 0 , 


T = e^^ + (1-a^^) e 


22 


+ a. 


\ 1-Uj 


•12 
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APPENDIX G 


Computer Program for Calculating 
Lunar Tidal Stresses and Strains 


due to the Earth 



n o n a o 


LUMAR stress aHD strain from earth tides 

M HOUSTON 5/lH/7(!! 

this program calculates the tidal stresses and strains 
ON the HOON due to THE EARTH, 

IMPLICIT double precision (A-H^L,0"I) 

dimension JIEAR (li)/26) ,JDAT(it»1^)^E(3,i),Q(3f3) 

data JTEAR(0)/241S019/^ ( ODA I ( 1 , K ) , 1 , 1 2 ) / 

n?/3l^60> 91, 121, 152, 182, 213,244,27 4, 30b, 335/ 

COSQ(THETA) = DC0S{THETA) >*DC0S(THETA) 

SINQ(THETA) = DSIN(THETA)*DSIN(THETA) 

C0S2(THETA) = DCOS(2.0D0>fTHETA) 

SIN2(THETA) =: DSIN(2.0D0H'THSTA) 

GHOON = 1.67D2 

EMUXST = 3.84411D10 

RMOOK = (2.73D-1)+(6,37123D8) 


400 


401 


402 


C 

C 

C 

1 


2 


99 


type 400 

F0RMAT(/' lunar love NUMBERS, K AND L V) 

accept 91/HLOVE 
accept 91,LL0VE 
type 401 

FQRMAT(/' LUNAR COMPRESSIBILITY AND RIGIDITY, LAMBDA MU'/) 

accept 91, lambda 
accept 91,LMU 

type 402 

P0RMAT{/' direction cosines for stress direction '/) 
accept 91#L1 
accept 91»L2 
accept 91, L3 

CONST = ( 3 . 0D0 ) fRMOON* ( 6 , 67D-8 ) * { 5 . 9B3D27 ) 

CONST = CONST/ ( ( 4 , 0D0 ) * EMDlST*EHDlST*EMDlST*GMOON ) 

JULIAN calendar 4/2/70 
DO 1 J=1,26 

JYEAH(J) = JYEAH(J-I) + 1461 
DO 2 1=2,4 
DO 2 0=1,12 

JDAY(I,0) = JDAY(I-1,J) + 365 
JDAY{ 2, 1 ) = JDAY(2, 1 ) + 1 
JDAY(2,2) = jDAYf2,2) + 1 
TYPE 99 

F0RMAT(/' number of CASES (12)'/) 
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NCASE = 

ACCEPT ICASF 

999 TYPE 100 

100 FORlUl't/' DAX/MOKTH, year - :H2'/) 

NCASE = NCASE + 1 

accept 10 i,iday,imonth,iyear 

101 f0REAT(3I2) 

11 = iYEAH/4 

12 = -Ut-II+IYEAR+l 

JULDAI = JXEAH(II) + JDAT(l2,IMaNTH ) H- IDAY 
type 102,JULDAY 

102 rORMAT{/' tJULlAN CALENDAR DAY =',H5) 

type 103 

103 FORmAT(/' NUHtiER OF DAYS -1 FOR WHICH RESULTS DESIRED = V) 
accept 104,NDAY 

104 F0BMAT(I2) 

hday = NDAY+jULDAY 

type 90 

90 F0RHAT(//' lunar OBSERVATION POINT'//' LATITUDE = '/) 

accept 91,0SLAT 

91 format (D) 
type 92 

92 FOEMAT(/t LONOITUDE = '/) 

accept 91(.0BLUNG 

type 200 

200 PORHAT(/t OUTPUT FILE UNIT DESIGNATION (12) “ '/) 

accept 134(,XUNIT 

type 201 

201 FORMAT(/i output FILE NAME (Ab) = '/) 

accept 202,ZNAMB 

202 FORMaT(A5) 

CALL 0FILE(IUNIT,ZNAME) 

WHITE (IUNIT, 190)OBLaT, OBLONG 

190 FOEmaT(///, IXJSHLUNAR observation POINT COORDINATES^// 

1 7XBHLATIIUDE, / X9 HLONGITU DE , // , 5 X 1 PD 1 3 , 6 , 2X , D 1 3 , 6 ) 


HDBG = 1 .745329251994329D-2 

oblat = ohlat+rdeg 

OBLONG = 0BL0NG*RDEG 
julday = JULDAY - 1 
fi9 JULDAY = JULDAY + 1 
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c 

c 


1.UHAR EPHEHEHIS 4/;>0/7H 


T = DBLE{ELOA’i{JiiLDAX-24l5019) )/3.6b2SD4 

pj, = i, 14 1592U53589793D0 

Al = (9.21D1)*Pl/( ( 1.8D2)*(6.0D1 ) ) 

S = 2.7043659U2 + ( { ( 2 . 0D-6 ) #'X'+ i . 98 3 ) * T + 4 , 8 1 2 6 7 89057 DS ) * T 
P = 3.34329561^2 + ( { ( - 1 , 0D-5 ) *T^ 1 , 032D-2 ) *T + 4 . 06903403D3 ) *T 
H = 2.79696681^2 + (( 3 .00-4 ) *T'-3 . 600076892^4 ) 

RfllN = 2,9088820866572150-4 
S = S*R0E9 
P = P*RDE0 
H = H*'E0EG 


DRATIO p 1.0D15 + (5.45D-2)*0C0S(S-P) + C 2 . 9 70- 3 ) *DCOS { ( 2 , 0D0 ) * 
1(S^P)) + ( 1 ,0020-2) fDCOS{S-P-2.0D0* (H^P) ) 

2 + (8,250-3) *DCOS(2,0O0'f (S-H ) ) 

pS = 2.8 12208302 + ( { ( 3 , 00-6 ) f T + 4 . 50-4 ) *T+ 1 . 7 1 90200 ) *T 
pS = PS>t<SDEG 

51 = (3,7/02) *DSIN(S-P) + ( 1 . 3D 1 } +DsiN ( ( 2 , 000 ) * ( S-P ) ) 

H- {7.6D1)^DSIN(S-P-(2,0D0)>K(H-P) ) - ( 1 . 1 D 1 ) *DSIN ( H-PS ) 

2+ (4.0O1)*O3IN( (2,0D0)*(S-H) ) 

Si = 5 H< SUEMIH 


AN = 2,591832802 + { ( ( 2 . 0D- 6 ) *T (■ 2 , 08 0-3 ) *T- 1 . 9 3 4 1 4 20 1 0 3 ) *T 
an = AR*RDEG 

A = DSIN( AI) *0C0S{S1_AN) 

AMU = (1.03D-2)*DSIR((2,0D0)*(S1-AN) )#ROBG 


beta = (3.08D2)*dSIN(S-AN) + ( 1 . 70 1 ) fOSIN ( 2 , 0D0*S-P-AN ) 
1+ ( 1 .001) ♦DSIN{ i2.0D0)*(S-H)-S + AN) - ( 1 . 70 1 ) *DSlN ( P- AN ) 

beta = BETAkRMIN 

TANb = -DSIN(AI)*DS1N(S1-AN)/DCoS(AI) 

ELAT = (DATAH('I'ANB)-BBTA) 

ELONG = (SI-S-1-AMII + A + ELAT) 

TANB = -DSIN(AI)>kDSIN(S1-AN)/DC0S(AI) 

ELAT = (DATANITaNB)-BSTA) 

ELONG = (81-S+AHU+A*ELAT) 
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2fe54 

C 


K1 ■= iv'LAT/RDE^ 

VU = KLONG/RD^'-IG 

WRITE (lURlT, 30J \JULD<lir 

P0RMAT(/////1Jt4feHBARTH SELENOGRJVPHIC COORDIHATES FOR JULIAN DAY, 

1 114 ) 

WRITE (IUN1T,^W4)E1,E2,DRATI0 

FORM at ( //»VX8HLATITUDB,7X9H longitude, fiXBHPARABL AX, 

! //,5X, 1PD13.6,2(2X,D13.6) ) 

STRESS AitU strain CALCULATIONS 7/21/70 

DRATIO = UHATIO*DRATIO>kI)RATXO 
P - OiSLAT 
1) - ELAT 

L = ELONG - OBLONG 

E(3,3) - COSQ(P)*COSt!(D)>fCOS2(L)+SlN2(P)*SlN2(D)#DCOS(L) 

1+ ( (3,fe)D0) +SINQ(P)-1 .000) * (SINQ(D )- 1 .0D0/3.0D01 

CUBIC = {4.0D^5>t'HLOVE«6,f'5D0*LLOVE) *C0NST.(.E(3,3) *DRATIC 
E(3,3) - (2.aD0)x.HLOVE*CONST*E(.'i,3)*DUATIO 

E(2,2) = 2.0Dl'*LlOVE>!‘{-3IN2(P)S‘CO32{p)*(COSQCD)*COS2(L) 

1 >2,0O0*COS2(D)’«DCOS(L) } +3 . 000* SIN 2 ( P ) * SlNQ ( D ) -SIN 2 ( P ) ) 

2 -f-HLOVE*(COSQ(p) *C0SQ(D)*C0S2{L)+SIN2( p) *SIN2{D) *DCOS (L) 

3 ^-(3.0O0*SINQ(p;-1.01J0)*(SINQ(D)-1.0D0/J.0D0) ) 

E(2,2) = C0NST*E(2,2)*DRAT10 

E{1,1) = -{LLUVE/COSQ(P} )*(C0SQiP)*C0SQ(D)*CDS2(L) 

1 4-SIN2( P) *SIN2 ( U) #DCOS(L) ) - ( LL0VE*DSIN ( P ) /DCOS (P ) ) * 

2 (SIN2{?) *COSQ(U) *cOS2(L)+2.0O0*CnS2(P) * S I N 2 ( D ) * DCOS ( L ) 

3 +3.0D0+SIN2(P) *SIKQ(D)-SIN2(P) ) + HLOVE* ( COSQ ( P 1 *COSQ ( D ) *C0S2 ( L ) 

4 *SIN2(P) *SXN2(0) *DCOS(L)4.(3.!0D(/)*rlNQ(P)- l .000)* (SINQ(D) 

3 -1 .000/3,000 ) ) 

E{1,1) = CONS'i'*E( 1, 1 ) *DHATI0 

E(1,2) = 4,gD«'J*(PSIN(P)*COSQ{D)*SlN2(L)-DCOS(P)*SXN2{D)*DSlN(M) 
I ♦>SIN2(0)*DSIN(L)/DCOS(P) 

E(1,2) = LLOVt;*CONST*B( 1,2 )*ORA'j:IO 

E(1,3) = 2 .00 J*HLOVE# (-DCOS ( p ) *C0SQ ( D ) *SIN2 CL) 

1 -5,0D-1*DS1N(P)*SIK2{D)*DSIN(L) 

2 H-COSQ(P) *COSQ(U)#COS2(L)+SIN2(P)«SIN2(D) *DCOS(L) 

3 '■(3.0D0*SINQ(P)-1 ,p;D0) *(SINQ(C .000/3,000) ) 

4 V LL0VB*(DC0S(P)*C0SQ(D)>rSlN2(L) 

3 - 2.0D0*DSXN(p ) * sIN2(D)*USIN(L) ) 

EC 1,3) = CONS'i'*EC 1/3)*Dr'ATX0 

r(2,J) =. CHLOVE + LLOVE)*C-SIN2 (P)*COsqCD)*COS2CL) 
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1 <-2,0D0*COS2(P)*SIN2(D)«kDCOS{L) 

2 +3.0D0*SIN2(P) *SINQ(D)-SIM2(P) ) 

E(2,J) = CONSl'*Ef2,3}*DI!ATlO 

WHITE (IUNIT,403)E( 1,1),E(1,2)»K(2#2),E(1,3),E(2,3).E{3,3) 
F0RHAT(///, U I9HTIDAL STRAIN TENSOR ,//^ 

1 10X3HKXX,/^5X iPD'i3 .6,//, 10X3HEXY, 1 2X3HEY Y , / , 5X 

2 2(D13.6,2X),//, 10X3HEXZ, 12X3HEY:-, 1 2X3H EZZ , / , 5X , 3 ( D 1 3 . 6 , 2 X ) ) 

TIDRL STREiiS CALCULATION 
DO 450 I='i,3 

Q(1,I) = LAHBOA*E(I,l ) + 2,0D0*LMU*E(I^I) 

0(1,2) = 2.0d0*LMU*E( 1,2) 

0(2,3) = 2.0D0*LMUt‘E(2,3) 

0( 1,3) “ 2,0d0*LMU*E( 1,3) 

WRITE (IUNIT,404)O(1,1),Q(1,2),Q(2,2),Q(1,3),Q(2,3),g(3,3) 
F0RMAT(///, 1X19HTIDAL STRESS TENSOR ,//, 

1 10X3HTXX,/,5X1PD13.6,//, 10XJHTX1, 1 2X3 HT Y Y , / , 5X 

2 2(D13.6,2X ) ,//, 10X3KTXZ, 12X3HTYZ, 1 2X3HTZZ , / , 5X 3 ( D 1 3 . 6 , 2X ) ) 

WHITE ( IUNIT,405)CUBIC 

FORMaT(//,' cubic dilation = S'lFDli.b) 

TD = L1*L1»Q{ 1,1 )+L2*L2 + Q{2,2)+L3^L3*g(3,3)+L 1*L2*0( 1,2) 

NOTE THAT BSCaUSE oF LIMITING CONDITIONS, NEAR THE SURFACE 
OF THE MOON, T(1,3) = T(2,3) 0 


WRITE (IUNIT,406)TD,l1,L2,L3 

F0RMAT{//, 1X6i1STRESS , 9X 1 7HDIRECT10N COSINES , /2 ( 1 PD 1 3 . 6 , 2X ) , 
1 2(/15X,Dl3.6) ) 

IF (JULDAX-NDAY) 69,70,70 

CONTINUE 

end file IUNIT 

IF (NCASE-ICaSE) 999,71,71 

CONTINUE 

type 302 

FORMaT(//' FINIS') 


STOP 

end 
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